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Recent developments in quantum physics make heavy use of so- 
called "quantum trajectories." Mathematically, this theory gives rise 
to "stochastic Schrodinger equations," that is, perturbation of 
Schrodinger-type equations under the form of stochastic differential 
equations. But such equations are in general not of the usual type as 
considered in the literature. They pose a serious problem in terms of 
justifying the existence and uniqueness of a solution, justifying the 
physical pertinence of the equations. In this article we concentrate on 
a particular case: the diffusive case, for a two-level system. We prove 
existence and uniqueness of the associated stochastic Schrodinger 
equation. We physically justify the equations by proving that they are 
a continuous-time limit of a concrete physical procedure for obtaining 
a quantum trajectory. 

1. Introduction. Belavkin equations (also called stochastic Schrodinger 
equations) are classical stochastic differential equations describing the evo- 
lution of an open quantum system undergoing a continuous quantum mea- 
surement. The solutions of such equations are called quantum trajectories 
and describe the time evolution of the state of the system. The random na- 
ture of the result of quantum measurement is at the origin of the stochastic 
character of the evolution. 

The first rigorous description of a state undergoing a continuous measure- 
ment is due to Davies in [4]. It describes in quantum optics, the behavior of 
an atom from which we observe the photon emission. This is the so-called 
"resonance fluorescence" experiment (see [6] and [2]). 

In the literature, essentially two kinds of Belavkin equations are consid- 
ered: they are driven either by a Brownian motion or by a counting process. 
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But the kind of equations which are obtained this way are of nonusual type 
compared to the usual theory of stochastic differential equations. In partic- 
ular there is no reference in physics nor in mathematics, where the existence 
and the uniqueness of the solution of such equations are discussed. Further- 
more, the physical justification of the apparition of these equations requires 
in general a quite heavy mathematical framework (Von Neumann algebra, 
conditional expectation, filtering, etc.). The high technology of such tools 
contrasts with the simplicity and the intuition of the physical model. 

An approach to such equations, which is physically very intuitive, is the 
one of repeated quantum interactions. The setup is the following. The con- 
tinuous measurement model is obtained as a limit of discrete models. This 
discrete model is a naive approach to the interaction of a simple system 
interacting with a field. The field is represented as a chain of independent 
copies of small pieces of environment. The simple system interacts, for a 
time interval h, with one piece of the environment. After that interaction 
an observable of the piece of environment is measured. The random result 
of the measurement induces a random new state for the small system. The 
small system then interacts again with another piece of the environment for 
a time interval h. A measurement of the same observable of this second copy 
is performed. And so on. 

This experiment gives rise to a discrete evolution of the state of the small 
system, which is a Markov chain. The continuous-time limit (h — ► 0) of this 
evolution should give rise to the quantum trajectories. 

Repeated quantum interactions have been considered by Attal and Pau- 
trat in [1] and by Gough in [5]. The continuous limit of repeated quantum 
interactions is rigorously shown to converge to a quantum stochastic differ- 
ential equation in [1]. The setup of measuring an observable of the chain 
after each interaction is considered in [5], but the continuous limit, the ex- 
istence and the uniqueness of the solutions are not all treated rigorously in 
this reference. 

The aim of this article is to study the diffusive Belavkin equation, to 
show existence and uniqueness of the solution, to show its approximation 
by repeated quantum interaction models. The same results for the equation 
concerning the counting process are developed in another article [11]. 

This article is structured as follows: In Section 1, we present the math- 
ematical model of repeated quantum interactions with measurement. We 
define discrete-time quantum trajectories and focus on their probabilistic 
properties. In particular, it is shown that these processes are classical Markov 
chains. Finally we deal with the model of a two-level atom in interaction with 
a spin chain and we describe the discrete stochastic evolution equations in 
this setting. 

Section 2 is devoted to the continuous model. We present the two different 
types of Belavkin equations whose solutions are continuous-time quantum 
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trajectories. We then prove existence and uniqueness of solutions in the 
diffusive case. 

The link between discrete and continuous models is provided in Section 
3. It is shown that particular discrete quantum trajectories (for a two-level 
model) satisfy stochastic equations which are discrete-time diffusive equa- 
tions. We use result of weak convergence of stochastic integrals in order to 
prove that solutions of diffusive Belavkin stochastic equations are obtained 
as a limit of discrete trajectories. 

2. Discrete quantum trajectories. We make here precise the mathemat- 
ical framework to describe the model of discrete quantum trajectories. 

2.1. Repeated quantum measurements. The physical setup is the one of 
a small quantum system, represented by a Hilbert space Tio, coupled to a 
field modeled by an infinite chain of identical independent quantum systems. 
Each piece of the field is represented by a Hilbert space Ti. Each space is 
equipped with a positive, trace-class operator with trace 1. This operator 
is called a "state" or "a density matrix." In this section, we present the 
random character of repeated measurements. 

The discrete model of interaction is called "quantum repeated interac- 
tions." Each copy Ti of the environment interacts with Tio, one after the 
other, during a time interval of length h. Information on the evolution of 
the small system is obtained by performing a measurement of Ti after each 
interaction. 

For the first interaction, the compound system is described by the tensor 
product Tio ® Ti and the interaction is characterized by a total Hamiltonian 
H tot which is a self-adjoint operator on Tio (g) Ti. Its general form is 

(2.1) H tot = H ®I + I®H + H hlt , 

where Hq and H are the free Hamiltonians of each system and H{ nt is the 
interaction Hamiltonian. The operator 

jj — gihHtot 

describes the first interaction as follows. In the Schrodinger picture, if p 
denotes any state on the tensor product, the evolution is given by 

p^UpU*. 

After this first interaction, a second copy of Ti interacts with Tio in the same 
way. And so on. 

As the field is supposed to be an infinite chain, the whole sequence of 
successive interactions is described by the state space 

(2.2) T = Ho®<g)H k , 

k>l 
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where Tik designs the fcth copy of Ti. The countable tensor product <S>k>i %k 
is defined as follows. We consider that Ti is finite dimensional and that 
{Xq,X\, . . . ,X n } is a fixed orthonormal basis of Ti. The orthogonal projec- 
tor onto CXq is denoted by |Xo)(Xo| (this is the braket notation in mathe- 
matical physics; see the remark below). This is the ground state (or vacuum 
state) of Ti. The tensor product is taken with respect to Xq (for details, see 
[1]), that is, we define an orthonormal basis of (8)/c>i ^fc with respect to this 
vector. It is described as follows. 

Let V be the set of finite subset A of the form A = {(ni, . . . , (n^, ik)} of 
N* x {1, . . . ,n} such that the re^'s are two-by-two disjoint. The orthonormal 
basis of ®fc>iWfc with respect to Xq is the family 

{X A ,AeV}, 

where for A = {(ni,ii), . . . , we define Xa as the vector 

Xq <g) ■ ■ ■ <g) X h (g> Xq <g) ■ ■ ■ <g) Xq <g> X h ® ■ ■ ■ , 

of ®fc>i7^fej where Xj. appears in the copy number nj of Ti. The infinite 
tensor product allows us to work in a single space but the structure of Hilbert 
space does not appear explicitly in the rest of the paper. 



Remark. A vector Y in a Hilbert space Ti is represented by the appli- 
cation \Y) from C to TC which acts with the following way: |^)(A) = |AY). 
The linear form on Ti is represented by the operators {Z\ which act on the 
vector \Y) by (Z||Y) = (Z,Y), where (•, •) denotes the scalar product of Ti. 

The unitary evolution describing the fcth interaction is given by the oper- 
ator Uk which acts as U on TiQ®Tik, whereas it acts as the identity operator 
on the other copies of Ti. If p is a state on T, the effect of the kth interaction 
is 

p^UkpUt. 

Hence the result of the k first interactions is described by the operator Vk 
on /3(r) defined by the recursive formula 

/ 9 / Vk+i = Uk+iVk, 

[2 - 6} \Vq = I, 

and the evolution of states is then given, in the Schrodinger picture, by 
(2.4) P^V kP V k \ 



We present now the indirect measurement principle. The idea is to perform 
a measurement of an observable of the field after each interaction. 
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A measurement of an observable of 7ik is modeled as follows. Let A be any 
observable on TC, with spectral decomposition A = Ysj=i ^jPj- We consider 
its natural ampliation on T: 

fc-i 

(2.5) A k :=(g)I®A® /. 

3=0 j>k+l 

The result of the measurement of A k is random; the accessible data are its 
eigenvalues. If p denotes the reference state of T, the observation of Aj is 
obtained with probability 

P[to observe Xj] = Tr\pl*], j & {1, . . . ,p}, 

where P k is the ampliation of Pi in the same way as (2.5). If we have observed 
the eigenvalue Xj, the "projection postulate" called "wave packet reduction" 
imposes that the state after measurement is 

PfpPf 
Pj ~ Tr[pP/] ' 

Remark. This corresponds to the new reference state depending on the 
result of the observation. Another measurement of the observable A k (with 
respect to this new state) would give P[to observe Xj] = 1 (for P%Pj = if 
i j). This means that only one measurement after each interaction gives 
significant information. This justifies the principle of repeated interactions. 

The repeated quantum measurements are the combination of the previous 
description and the successive interactions (2.4). After each interaction, the 
measurement procedure involves a random modification of the system. It de- 
fines namely a sequence of random states which is called "discrete quantum 
trajectory." 

The initial state on T is chosen to be 

3>1 

where p is any state on Hq and each = is a fixed state on TC. We denote 
by pk the new state after k interactions, that is, 

Pk = VkpVt- 

The probability space describing the experience of repeated measurements 
is f^* where O = {1, . . . ,p}. The integers i correspond to the indexes of the 
eigenvalues of A. We endow N with the cylinder u-algebra generated by 
the sets: 

Aj lv ..,j fc ={w£Sl N /^i =h,...,u}k = ik}- 
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Note that for all j, the unitary operator Uj commutes with all P k , for k < j. 
For any set . . . ,i k }, we can define the following nonnormalized state: 

. . . , i k ) = (I ® P h ® • • ■ g> P ik ® I ■ ■ -)nk{I O -P*! ® • • ■ ® P< fc ® / • ■ •) 

= { p i k --- p l)Vk{Pl 1 ---Pi k )- 

It is the nonnormalized state which corresponds to the successive obser- 
vation of the eigenvalues , . . . , \i k during the k first measurements. The 
probability to observe these eigenvalues is 

P[Ai lv ..,jJ =P[to observe (X h , . . . , X ik )\ = Tr[/i(ii, . . . ,i k )]. 

This way, we define a probability measure on the cylinder sets of f2 N * which 
satisfies the Kolmogorov Consistency Criterion. Hence it defines a unique 
probability measure on N . The discrete quantum trajectory on T is then 
given by the following random sequence of states: 

p k --n n * — >B(r), 

- / \ Ji{u-i, . . . ,u k ) 

UJ I > Pk{0Ji,...,U k ) = — -— -r. 

Tr[fi(ui,...,u k )\ 

From the construction and the remarks above, the following is immediate. 

Proposition 2.1. Let (p k ) be the above random sequence of states; we 
have for all u) G J7 N : 

Pk+l[UJ> Tr{p k (u)U* +1 P^U k+1 ] ' 

The following theorem is an easy consequence of the previous proposition. 

Theorem 2.1. The discrete quantum trajectory (p n ) n is a Markov chain, 
with values on the set of states of Ti.Q^ i>1 TCi. It is described as follows: 

P[Pn+l = P/Pn = n , ■ ■ ■ , po = 9 ] = P[p n +1 = p/Pn= On]- 

If Pn — On, the random state p n +i takes one of the values: 

pp + \u n+1 (o n ^p)Un + i)pr +1 



Tr[(U n+1 (0 n ®(3)UZ +1 )Pl 



n+l] 



1,...,P, 



with probability Tr[(U n+1 (O n ® /3)U* +1 ) P? +1 ] 



In general, one is more interested in the reduced state on the small system 
Tio only. This state is given by taking a partial trace on Hq. Let us recall 
what partial trace is. If Ti is any Hilbert space, we denote by Tr-^fVF] the 
trace of a trace-class operator W on 7i. 
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Definition 2.1. Let TC and /C be two Hilbert spaces. If a is a state on 
a tensor product TC <8> /C, then there exists a unique state 77 on TC which is 
characterized by the property 

Tr H [r j X]=Tr 7mlc [a(X ® I)} 

for all X 6 B(TC). This unique state 77 is called the partial trace of a on TC 
with respect to K,. 

Let Eo(a) denote the partial trace on TCq with respect to <S>fc>i'^fc °f anv 
state a on T. We define a random sequence of states on TCq as follows. For 
all oj in f2 N , define the discrete quantum trajectory on TCq 

(2.6) p n (cj) = E [p n (u/)]. 

An immediate consequence of Theorem 2.1 is the following result. 

Theorem 2.2. The quantum trajectory {p n ) n defined by formula (2.6) 
is a Markov chain with values in the set of states on TCq. If p n = Xn, then 
p n +i takes one of the values: 

F \ (I®Pi)U{Xn®P)U*{I®Pi) } 

with probability Tr[U(x n ® P)U*(I <g> Pi)} . 
Remark 1. Let us stress that 

{I®Pi)U{Xn®P)U*(I®Pj) 

Tr[U(xn®P)U*(I®Pi)\ 

is a state on TCq (8) TC. In this situation, the notation Eo denotes the partial 
trace on TCq with respect to TC. The infinite tensor product T is just needed 
to have a clear description of the repeated interactions and the probability 
space fi N . 

The next section is devoted to the particular case of a two-level atom in 
contact with a photon stream. Because of physical considerations, this case is 
often the central case in the literature concerning continuous measurement. 

2.2. A two-level atom. The Hilbert spaces describing the physical situa- 
tion are now TCq = TC = C 2 . 

In this section, we establish a discrete quantum evolution equation for 
(Pn) which is a discrete approximation of the Belavkin equation. 

The main goal of this section is to obtain a formula of the following form: 



(2.7) 



Pk+l — f(Pk,X k+ i), 
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u 



where (Xk)k is a sequence of random variables. Such a formula is obtained 
from the description of Theorem 2.2 and the computation of the partial 
trace operation. 

The state pk can be considered as an initial state (according to the Markov 
property of Theorem 2.2). Thus we can consider a single interaction with a 
system (H,P) [actually this is the (k + l)st copy]. We consider an observ- 
able of the form A = Ao-Po + A1P1 and the unitary operator describing the 
interaction is a unitary 4x4 unitary matrix. 

In order to compute the state given by the projection postulate and the 
partial trace, we choose a suitable basis. If (Xq = Q,Xi = X) is an or- 
thonormal basis of C 2 , for the space TLq <g> 7i we consider the following basis: 
f2(g)J7,X®f2,f2(g>X, X <S> X. In this basis, the unitary operator U can be 
written by blocks in the following way: 

^00 Lqi 

where each Ljj is an operator on TLq. For j3 we choose 

= \n)(n\. 

As a consequence, the state after the interaction is 

(, 8) 

Thanks to the description of Theorem 2.2, we define the two possible 
nonnormalized states 

(2.9) £ (p k ) = E [J® Po(p k +i)I ® Po], 

(2.10) d(p k ) = E [I® Pi(p k+1 )I O Pi]. 

These are operators on 7^o; the nonnormalized state £o(Pfc) appears with 
probability pfc+i = Tr[£ (Pifc)] and£i(/9fc) with probability g fc+ i = Tr[£i(/9 fc )]. 
Let us define the random variable Vk+i on {0, 1}: 

f ffc+i (0) = 0, with probability , 
\ ffc+i(l) = 1, with probability q k+1 . 

With these notation, the discrete quantum trajectory can be described as 
follows. For all u> G , we have 

(2.11) p fc +i(w) = -— (1 - ffe+i(w)) + — — 

Pfc+iM 

In order to obtain the final discrete quantum evolution equation, we con- 
sider the centered and normalized random variable 
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We define the associated filtration on {0, 1} N : 

T k = a(Xi,i < k). 

By construction, we have E[Xk + i/Fk] = and E[X% +1 /Fk] = 1. Thus we can 
write the discrete evolution equation for the quantum trajectory in terms of 
the random variables (X).): 



(2.12) Pk+l = C {p k ) + C 1 {p k ) + 



Pk+i 



Qk+i 



X 



fe+i- 



The above equation can be considered in a general way and the unique 
solution starting from po is the discrete quantum trajectory described in 
Theorem 2.2. Let us stress that this sequence depends on the length of 
time of interaction. This dependence will allow us to prove a continuous- 
time approximation result in Section 3. For the moment, the next section is 
devoted to describing continuous-time quantum trajectories. 

3. Belavkin equations. As was announced in the Introduction, it is com- 
monly assumed that the evolution of a system undergoing a continuous 
measurement is described by stochastic differential equations. A model of 
interaction can be provided to describe an atom in contact with a continuous 
field. In this setting, the description of the principle of indirect measurement 
needs highly technical tools in order to obtain rigorous statements. Such the- 
ories are not the purpose of this article. We just give the physical setup in 
order to introduce the Belavkin stochastic differential equations. 

Consider a two-level system, described by C 2 , in interaction with an envi- 
ronment (classically described by a Fock space). The time evolution is given 
by a unitary process (Ut) which satisfies a quantum stochastic differential 
equation (cf. [10]). Without measurement the evolution of the small system 
is given by a norm-continuous semigroup {Tt}t>o- The Linblad generator of 
(Tj) is denoted by L and we have the "Master equation": 

^ = L( Pt ) = -i[H , Pt ] - \{CC\ Pt } + C Pt C\ 

where C is any operator on C 2 and Hq is the Hamiltonian of the atom. 

In the theory of time-continuous measurement L is decomposed as the 
sum of C + J where J represents the instantaneous state change taking 
place when detecting a photon, and C describes the smooth state variation 
in between these instants. These operators are defined by 

£(p) = -i[H ,p]-±{CC\p}, 

J(p) = CpC\ 

Thanks to the work of Davies in [4], two types of stochastic differential 
equations can be derived. The solutions of these equations are then called 
"continuous-time quantum trajectories": 
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1. The "diffusive equation" (homodyne detection experiment) is given by 
(3.1) d Pt = L(p t ) dt + [p t C* + C Pt - Tr( ft (C + C*))p t ] dW t , 

where Wt describes a one-dimensional Brownian motion. 

2. The "jump equation" (resonance fluorescence experiment) is 

(dN t -Tr[J(p t )]dt), 

where N t is a counting process with stochastic intensity J Q Tr[J~(p s )] ds. 

A physical justification of (3.1) as limit of discrete quantum trajectories 
is given in Section 3. For the moment, we shall focus on the general problem 
of existence and uniqueness of a solution of (3.1). The jump equation and 
all the convergence theorems referring to this case are treated in detail in 
[11] with different techniques. 

3.1. Existence and uniqueness. Let po be any state; we aim to show 
existence and uniqueness for the stochastic differential equation 

(3.3) Pt = Po+ f L(p s )ds+ [ t [p s C* + Cp s -Tr(p s (C + C*))p s }dW s . 

Jo Jo 

Classical theorems concerning existence and uniqueness for SDE cannot be 
applied directly here for the coefficients of (3.3) are not Lipschitz. Further- 
more, even if there exists a solution, one must show that the solution takes 
values in the set of states. Actually this property and the questions of exis- 
tence and uniqueness are linked. 

Concerning the property of being valued in the set of states, an important 
feature of the differential (3.3) is that it preserves the property to be a pure 
state (in quantum theory, a pure state is a one-dimensional projector). This 
idea is expressed in the following proposition. 

Proposition 3.1. Let (Wt) be a standard Brownian motion on (£l,J-,J-t, 
P) and let \ipo) be any norm-1 vector in C 2 . Let vt = ^(ifH, (C + C*)tpt) where 
C is any operator on C 2 . 

If the following stochastic equation 

(3.4) d\iln) = (C- v t I)\^ t ) dW t + (-m - l(C*C - 2u t C + vp))\i> t ) dt 

admits a solution (\ipt))> then almost surely we have \\ipt\\ =1 for all t. 

Furthermore the process (\iftt) {ipt\) takes values in the set of pure states 
and it is a solution of the diffusive Belavkin equation (3.3). 



(3.2) 



dpt = L(p t )dt + 



JiPt) 
TrL7(pt)l 



Pt 
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Proof. Let \ipo) be any vector in C 2 and let (\ipt)) be a solution of (3.4). 
Let us prove that \\ipt\\ 2 = 1- Using the Ito formulas and the fact that H is 
self-adjoint, a straightforward computation shows that 

d\\if)t\\ 2 = d(ipt,ipt) = (dipt,ipt) + (A, dipt) + (dipt, dipt) 
= {{C-v t I)Mt)dW t 

+ {(-iH - ±(C*C - 2v t C + v 2 I))^ t ) dt 
+ (iPt, (C - v t I)ip t ) dW t 
+ (rp t , (-iH - ±(C*C - 2v t C + u 2 t I))^ t ) dt 
+ {{C-u t I)^ u (C-u t I)^ t )dt 
= {2u t -2v t {^ t ))dW t . 
If HV'oll 2 = 1) this implies that almost surely 

W| 2 = W| 2 = l 

for all t > 0. Define the process pt = \ipt)(ipt\- It is valued in the set of pure 
states. As ||^|| = 1, we have for all y 6 C 2 

pt\y) = (ip t ,yMt)- 

Hence we can compute dpt\y) by the Ito formula: 

dpt\y) = (d?Pt,y)\il>t) + (^t,y)d\ip t ) + (difj t ,y}d\ip t } 

= ((C-u t )Tp t ,y)\i> t )dW t 

+ ((-iH - \{C*C - 2v t C + v 2 ))^ u y)\^ t ) dt 

+ (tPt,y)(C-u t )\^ t )dW t 

+ (ip t ,y)(-iH - \{C*C - 2u t C + i%))\*Ih) dt 

+ ((C-v t )^ t ,y)(C-v t )\^ t )dt. 

Let us show that this corresponds to (3.3). It is clear that vt = ^(ipt, (C + 
C*)ipt) corresponds to the term hTr[\ipt)(iJ)t\(C + C*)]. As a consequence 
the term in front of the Brownian motion becomes 

<(C - vt)^y)W) + &t,y)(c - v t )\rl>t) 

= (C'|^)(^| + |^)<^|C*-Tt[|^)(^|(C + C*)]|^)(^|)|y). 
A similar computation shows that the term in front of dt is 

m t )m\y)- 

Hence we recover the expression of Belavkin equation (3.3) and the propo- 
sition is proved. □ 
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As a consequence, we can express an existence and uniqueness theorem 
for (3.3). In what follows, we use the notion of "wave function." A wave 
function is a norm-1 vector which defines a pure state. 

Theorem 3.1. Let (fl,J-,J-t,P) be a probability space which supports a 
standard Brownian motion (Wt) and let po be any state on C 2 . 
The stochastic differential equation 

pt = p + f t L(p s )ds+ f'lpsC* + Cp s -Tr[( Ps (C + C*))p 8 ]]dW 8 
Jo Jo 

admits a unique solution (pt)- The solution takes values in the set of states 
and is defined for all t>0. 

Furthermore, if the initial condition is a pure state, the solution takes 
values in the set of pure states. The corresponding stochastic differential 
equation for a wave function is then given by 

d\i> t ) = {C- ut)\ih) dW t + {-iH - \{C*C - 2v t C + i/ 2 )M) dt 
where u t = l{ij H ,{C + C*)if>t). 

Proof. As the coefficients of (3.3) are not Lipschitz, we cannot apply 
directly the usual existence and uniqueness theorems for SDE. However, 
the coefficients are C°°, hence locally Lipschitz. We can use a truncature 
method. Equation (3.3) is of the following form: 

(3.5) dpt = L(p t )dt + e(p t )dWt 

where 9 is C°° and 0(A) = AC* + CA- Tr[A(C + C*)]A. Let k>l be an 
integer; we define the truncation function <p k from R to K defined by 

{—k, if x < —k, 
x, if —k < x < k, 
k, if —k < x < k. 

For a matrix A = (flij), we define by extension (f>k{A) = </?fc(Re(ajj)) + 
i(pk(Jxa(aij)) . Thus the function 9o^ is Lipschitz. Now we consider the 
truncated equation: 

dp k)t = L o <pk(pk,t) dt + @o <Pk{Pk,t) dW t . 

The Cauchy-Lipschitz theorem concerning stochastic differential equations 
can be applied; there exists a unique solution t \— > pkj defined for all t. 
Besides the solution is continuous in time. 
Define the random stopping times 

T k = M{t,3(ij)/\Re(aij(p k>t ))\ = k or \ Im(a,ij(p k j))\ =k}. 
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As po is a state, we have \\po\\ < 1- Thanks to continuity, if k is chosen large 
enough, we have Tk > and for all t < Tk we have <fk(Pk,t) = Pk,t- Thus 
1 1 — ► pfc,t is the unique solution of (3.3) (without truncation) on [0,Tfc]. The 
usual method for solving an equation with non-Lipschitz coefficients is to 
put T = linifc Tk and to show that T = oo. 

In addition to the proof of existence of a solution, we must prove that the 
process is valued in the set of states. If v is any state, we have ||z/|| < 1, so 
IKu)| — 1- Hence if we prove that on [0,T2] the process (p2,t) is valued on a 
set of states, this would prove that T2 = 00 a.s. and we would have proved 
that there exists a unique solution valued in the set of states. Let us prove 
this fact. 

In the proof of the existence and uniqueness of a solution in the case of 
Cauchy-Lipschitz coefficients, the solution is obtained as the limit of the 
sequence 

p n +l(t) =Pn(0) + / Lo(p k {pn{s))ds+ I 6 o (p k (p n (s)) dW s , 
Jo Jo 



(3.6) 

Po(t) = P- 

With our definition of and L, if po is a state, it is clear that this sequence 
is self-adjoint with trace 1. These conditions are closed and at the limit the 
process is self-adjoint with trace 1. But the condition of positivity does not 
follow from such arguments. We shall prove it by other means. 
Consider the random time 

(3.7) T° = inf{i < T 2 /3X £ C 2 /(X, p 2jt X) = 0}. 

We have (X,pqX) > for all X, so by continuity we have (X,p2jX) > on 
[0, T°] which implies that p2 : t is a state for all t<T . 

If T° = T2 a.s., the result is proved. Otherwise, if we have T° < T2, 
then by continuity there exists X such that {X, p 2 T°X) = and for all 
Y (Y,P2,T°Y) > 0. This implies that p 2 t° is a pure state because we work 
in dimension 2. Let us denote by ip^o a vector of norm 1 such that P2,t° = 
I ) (V't I • Consider the equation 

d\i> t ) = {C- v t )\fo) dW t + (-W - i(C*C - 2u t C + v 2 t ))\^ t ) dt 

with ip T o as initial condition. The problem of existence and uniqueness for 
this equation is solved by a truncation method also. The fact that, if we have 
a solution, it is of norm 1 shows that the solution obtained by truncation 
(defined for all t) is actually the solution of (3.4). Proposition 2.1 and the 
uniqueness of p2,t on P^ ,^] show that the solution of (3.4) which satisfies 

= hM + T (C - v s )\iJj s ) dW s + (-iH - \{C*C - 2u s C + u 2 s ))\^s) ds 
Jt° 

defines a process (\ip t ) {ijj t \) equal to p2,t on [T°,T 2 ]. Hence the process ob- 
tained by truncation is valued on set of states and the theorem is proved. 
□ 
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3.2. Change of measure. At this stage, it must be said that the stochastic 
differential equation usually appearing in the literature is of the following 
form: 

(3.8) Pt = Po+ [ t L(p s )ds+ f t \p 8 C* + C Pa - r Tr\p a (C + C k )]]dW„ 

Jo Jo 

where 

(3.9) W t = W t - /V[/9 S (C + C*)]cfe. 

Jo 

Hence it seems to be rather different from (3.3). Actually the link between 
the two different equations is given by Girsanov's theorem (see [12]). 

Theorem 3.2. Let (Wt) be a standard Brownian motion on (Q,J-,J-t,P) 
and let H be a cadlag process. Let 

(3.10) X t = fll.da + Wt 

Jo 

and define a new probability by 

i-«<-f ^-ijM 

for some T > 0. Then under Q, the process {X-t) is a standard Brownian 
motion for <t <T. 

The link between the two equations (3.3) and (3.8) is then obvious. Let 
(pt) be the solution of (3.3) given by Theorem 3.1 on (ft,F,F t ,P). For some 
T > 0, define the probability measure Q by 

(3.11) g = exp ^ Tr[ Pt {C + C*)] dW s - \ Tr[ Pt (C + C*)] 2 <fe) . 

The above theorem claims that Wt is a standard Brownian motion under 
Q for < t < T. Hence (3.8) is the same equation as (3.3) up to a change of 
measure. In the following section, we show that the solution of (3.3) can be 
obtained as limit of discrete quantum trajectories. 

4. Approximation and convergence. 

4.1. The discrete approximation. In this section, we present a way to ob- 
tain the solution of the diffusive Belavkin equation (3.3) as a limit of concrete 
discrete quantum trajectories. Let us show that these discrete trajectories 
satisfy evolution equations which appear as approximations of stochastic 
differential equations. 
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In Section 1, we had the discrete quantum trajectories satisfying 
(4.1) p k+1 = C (p k ) + d{p k ) + 



Hence, we have 

Pk+i — Po 

(4.2) 



Pk+i 



i=0 
k 

Y,[£o(pi) + £\(pi) - Pi 

i=0 



£o(pfc) + 



Pk+1 



Qk+1 



£\{pk) 



x, 



fe+i- 




£o{pi) + 



l Pi+i 

Qi+i 



Let us introduce a time discretization. Consider a partition of [0,T] in 
subintervals of equal size 1/n. The time of interaction is supposed now to 
be h = l/n; the unitary operator depends then on the time interaction: 



U(n) 



Loo(n) L m (n) 
Lio(n) L\i{n) 



In [1], Attal and Pautrat have shown that the asymptotic of the coeffi- 
cients Lij(n) must be properly rescaled in order to obtain a nontrivial limit 
when n goes to infinity. Indeed they have shown that Vj nt i = ET^i (n) . . . Ui(n), 
which represents the discrete dynamic of quantum repeated interactions, 
converges to a process Vt solution of a quantum Langevin equation only 
if the coefficients Ly obey certain normalizations. When translated to our 
context, the results of [1] show that we should consider 



(4.3) 
(4.4) 



£oo (n) 
L w (n) 



n 



-%H a - -ca 



»i - 

n 



1 



C + o 



v 



Remember that the unitary operator is given by 



U(n) 



expl i-Htot 
n 



The corresponding Hamiltonian H tot which gives the expression U(n) is of 
the following form: 



Hi 



tot 



H ®I + I< 



1 




+ 



1 





1 



+ c* 



1 




+ o\ 



r 

.n, 
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where Hq is the Hamiltonian of the small system and C is any operator on 
C 2 . 

With the time discretization we then obtain 
(4.5) Pk+i{n) = C (n)(p k (n)) + £i(n)(p fc (n)) 



+ 



(4.6) 



Qk+i{n'j 



£o{n)(Pk(n)) + 



I Pk+i{n) 
Qk+i(n) 



Ci(n)(p k (ri)) 



\j Pk+i(ji) 
x X k +i(n). 

Remember that the sequence of random variables (X k (n)) is defined 



through the two probabilities: 
(4.7) 

By definition of (X k ) we have 



(pk+i = Tv[C (p k )}, 
\q k+1 = Tr[£i(p k )]. 



(4.8) X k (n)(i) 



Qk+i(n) 
V p k+ i(n) 

lpk+i(n) 



with probability p k+ \{n) if i = 0, 



with probability q k+ \{n) if i = 1. 



Each probability depends on the expression of d, which depends on the 
measured observable A = Ao-Po + AiPi. At the limit, the diffusive behavior of 
Pint] ( n ) i s depending on the comportment of (X k ), that is, on the observable: 

1. If the observable is of the form A = Ao(q q) + Ai( ®), we obtain the 
following asymptotic for the probabilities: 

Pk+i(n) = 1-- Tr[J( Pk (n))} +o(-), 



n 



1 



q k+ i(n) = - Tr[J (p k {n))} + o 
n 

The discrete equation then becomes 
p k +i(n) - Pk{n) 



1 



L(Pk(n)) + o 



+ -Pfc(n +o (1) 

LTr[J(p fe (n))] 

x ^ (n)p fe+ i (n)X fc+ i (n) . 
2. If the observable A is nondiagonal in the basis (£l,X), for the eigenpro- 



jectors, put Pq 



(POO P01) 



and P 



(loo 901)- W e have 

^10 <?ll' ' 



Pk+i = Poo H 7= 1V[p fc OoiC + pi C*)] + - Tr[p fc poo(C + C*)] + o 
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Qk+i =qoo + -^r= Tr[p k (q 01 C + qi C*)} + - Tr[p k q 00 (C + C*)] +o(- 
\ n n \n, 



The discrete equation here becomes 

(4.9) p k+ i - Pk = -L(pk) + o( -) 
n \n J 

+ [e w C> fc + e- i6 p k C* 



(4.10) 



TY[p fc (e* y C + e -«C*)]p* + o(l))^=X k+1 . 

in 



In this expression, the parameter 9 represents a kind of phase. It is real 
and depends on the coefficients of the eigenprojectors. If we put Cq = 
e ie C, the discrete (4.9) becomes 



-L(p k )+o(-) 

n \n J 



+ [CePk + PkCt - Tv[p k (C e + C* e )]p k + o(l)]-^X k+1 . 



Pk+i - Pk 



For each 6, we have similar expressions for discrete equations with differ- 
ent operators Cq. Let us stress that this parameter does not modify the 
expression of L. In the following, we deal with 9 = 0. 

In [11], it is shown that the case where A is diagonal gives rise to the 
jump-Belavkin equation at a continuous limit. In the following section, we 
show that the diffusive case is obtained as the limit of the discrete process 
which comes from the measurement of a nondiagonal observable. 

4.2. Convergence theorems. Before presenting the main theorem con- 
cerning the convergence of discrete quantum trajectories, we show a first 
result concerning the average of the processes. In order to avoid confusion 
between the discrete-time process (p k ) and the continuous-time process (pt) 
we write the discrete process (p k ) with the index on the top. 

Theorem 4.1. Let (Q,X) be any orthonormal basis ofC 2 . For all non- 
diagonal observable A, the deterministic function t E[pt n *l (n)] converges 
in L°°([0,T]), when n goes to infinity, to the function t \— >— > E[pt]. That is, 

sup ||E[pM (ra)] _ E[ps] || o. 

0<s<T n ^°° 

Furthermore the function 1 1— >— ► E [pt] is the solution of the Master equa- 
tion 



dvt = L(i>t) dt. 
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Proof. First of all, we show the second part of the theorem. We can 
consider the function 1 1 > > E[pt] because we have existence and uniqueness 
of the solution of (3.3). The process {pt) is integrable (because pt is a state 
for all t). It is obvious that this function takes also values in the set of states. 
As po is a deterministic state we must show 

(4.11) E[p t ]=p + [ L(E[p s ])ds. 

Jo 

We know that the process (pt) satisfies 

Pt = Po + [ L( Ps )ds+ [ [ Ps C* + Cp s -Tr(p s (C + C*)) Ps ]dW s . 
Jo Jo 

As (Wt) is a martingale and the process (pt) is predictable (for it is 
continuous), the properties of stochastic integral give 

E [ f [p s C* + C Ps - Tt{p s (C + C*))p s ] dW s ] = 0. 
.Jo 

Hence, we have, by linearity of L, 

E[p t ]=Po + fvlL{p s )]ds 
Jo 

= Po+ f L(E[p s ])ds. 
Jo 

We then have the integral form of the solution of the Master equation and 
the second part is proved. 

Let us show now the first part of the theorem. We shall now compare 
E[pt n *](n)] with E[pi]. As in the continuous case, the martingale argument 
is replaced by the fact that the process (Xk) is centered. Remember that we 
have 

B[X k+l ] = E[E[X k+1 /F k }] = 0. 

As a consequence, considering k = [nt] and taking expectation in the discrete 
equation, we have 

[nt]-l 

E[pW(n)]-p°= -L(B[p k (n)]) + o(- 



i=0 

This is a kind of Euler scheme and we can conclude by a discrete Gronwall 
lemma argument that we have 

S u P ||E[pM( n )]_E[p s ]||-H.O. 

0<s<t n ^°° U 

The average of the discrete process is then an approximation of the aver- 
age of pt- In [1], this result was shown in the case of repeated interactions 
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without measurement. It is a consequence of the asymptotic of the unitary- 
operator coefficients, so it justifies our choice of the coefficients. 

Concerning the convergence of the processes, the discrete process which 
is the candidate to converge to the diffusive quantum trajectory satisfies for 
k = [nt] 

[nt]-l / i \ [«*]-! 

pM-p° 



£ -L(/(n))+o - + J2 [ (/)+ O (l)]-=X m . 
i=o n KnJ i=0 v" 

Thanks to this equation we can define the processes: 

[nt] 

v n k=l 
[nt] 



V n (t) 



(4.12) 



Pn {t) = pW{n) 



n 

nt]. 



[nt]-l [nt]-l 

e n (t)= £ o - + E ^T^^- 

By observing that these four processes are piecewise constant, we can write 
the process (p n (t))t>o as a solution of a stochastic differential equation in 
the following way: 

(4.13) Pn (t) = Po + e n (t) + f L(p n (s-)) dV n (s) + f @(p n (s-)) dW n {s). 

Jo Jo 

We now use a theorem of Kurtz and Protter (cf. [8]) to prove the conver- 
gence. Let us first fix some notation. 

Recall that the square-bracket [X, X] is defined for a semimartingale by 
the formula 

t 



[X,X] t = X?-2 I X s _dX s . 
Jo 



We shall denote by 7t(V) the total variation of a finite- variation process V 
on the interval [0, t\. The theorem of Kurtz and Protter [8] (see [7]) that we 
use is the following. 

Theorem 4.2. Let W n be a martingale and let V n be a finite variation 
process. Consider the process X n defined by 

X n (t)= Po + e n (t)+ f L(X n {s-))dV n {s)+ f Q(X n {s-))dW n {s). 
JO Jo 
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Assume that for each t>0, 

supEp^,Wn] t ]<oo, 

n 

supB[T t (V n )} <oo, 

n 

and that (W n ,V n ,e n ) converges in distribution to (W,V, 0) where W is a 
standard Brownian motion and V{t) = t for all t. Suppose that X satisfies 

x t = x + f ' L{x s )ds+ [ e(x s )dw s 

Jo JO 

and that the solution of this stochastic differential equation is unique. Then 
X n converges in distribution to X . 



We wish to apply this theorem to the process (p n (t)) [(4-13)]. The first 
step is the convergence of W n in (4.12) to a Brownian motion. We need 
the following theorem (cf. [3, 9]) which is a generalization of the Donsker 
theorem. 



Theorem 4.3. Let (M n ) be a sequence of martingales. Suppose that 

= 



lim E 

n— >oo 



sup|M n (s) -M n {s- 

s<t 



and 



[M n ,M n ] 



t. 



Then M n converges in distribution to a standard Brownian motion. The 
conclusion is the same if we just have 

lim E[|[M n> M n ] t -t|] = 0. 



Back to our process W n , consider the filtration 

F? = o{Xi,i<[nt]). 



Proposition 4.1. Let (W n ,V n ,s n ) be the processes defined in (^.12). 
We have that (W n (t)) is a J 7 ™ -martingale. The process (W n ) converges to a 
standard Brownian motion W when n goes to infinity. Furthermore we have 

supE[[W n ,Wg t ]<oc. 

n 

Finally, we have the convergence in distribution for the processes (W n ,V n ,e n ) , 
when n goes to infinity, to (W, V, 0) . 
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Proof. Thanks to the definition of the random variable X^, we have 
[Xi+i/Tf] = which implies E[^i" 
t> s, we have the martingale property 



E[X i+1 /Jf] = which implies E[± e£L ]+ i V^T] = for t > s. Thus if 



E[W„(t)/.??] = W„(*)+E 



1 



[nt] 



)1 



i=[ns]+l 

By definition of [Y, Y] for a stochastic process we have 



W n (s). 



nt] 



[W n ,W n ] t = W n (t) 2 -2 [ W n (s-)dW n (s) = -f]Xf. 

Jo n . =1 



Thus we have 



[nt] [nt] 

B[[W n , W n ] t ] = -E E[A, 2 ] = - -£B^[Xf/a{X h l < i}]} 

Tt . 

1=1 



n -I- 



iS 1 = M. 



i=l 



n 



Hence we have 



supE[[^ n ,W n ]i] <t <oo. 



Let us prove the convergence of (W n ) to (W). According to Theorem 4.3 
we must prove that 

lim B[\[M n ,M n ] t -t\}=0. 

n— >oo 

Actually we prove an Li convergence: 

lim B[\[M n ,M n ] t -t\ 2 ]=0, 

n— >oo 

which implies the L\ convergence. In order to show this convergence, we use 
the following property: 



B[Xf]=B[B[Xf/a{X l ,l< l }]] = l, 



and if i < j , 



E[(Xf - 1)(A| - 1)] = E[E[(X 2 - 1)(A| - l)/a{X h l < j}]] 
= B[(Xf-l)]B[(X]-l)} 
= 0. 



22 C. PELLEGRINI 

This gives 



E 



K^]t- — 
n 



•I \ 0-i -. [nil 



^Ee[(4 - 1) 2 ] + ^Eep, 2 - i)(*f - 1)] 

i=l i<j 
1 [nt] 



n 2 



i=l 



Thanks to the fact that poo an d qoo are not equal to zero (because the 
observable A is not diagonal!) the terms E[(X 2 — l) 2 ] are bounded uniformly 
in i so we have 

H N 2 



lim E 

n— too 



0. 



[Wn.W^ 

n 

As M— tin L „ 

we have the desired convergence. Finally, the convergence 
in distribution of iW n ) and (V^) implies the convergence of (e n ) to 0. □ 

We can now express the final theorem. 

Theorem 4.4. Let (£l,X) be any orthonormal basis of C 2 and let A be 
any nondiagonal observable (in this basis). Let po be any initial state on C 2 . 
Let (p^(n)) be the discrete quantum trajectory obtained from the quantum 
repeated measurement principle with respect to A. The process 

(pH( n )) then 

satisfies 

[nt]-l [nt]-l 



pW(n)=p + £ -L(p k (n)) + o(-)+ E [G(p fe )+o(l)]— X i+1 . 

i=o n Vny i=0 v re 

Let (pt) be the solution of the diffusive Belavkin equation (3.3) which satisfies 

Pt = po+ f L( Ps )ds+ f e( Ps )dW s . 
Jo Jo 

Then we have the convergence in distribution 



Proof. It is a simple compilation of Theorems 4.2, Proposition 4.1 and 
existence and uniqueness of Theorem 3.1. □ 
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